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ABSTRACT:  This  work  is  part  of  an  ONR  program  to 
imbed  DSP  in  next  generation  GPS  receivers  to  mitigate 
GPS  vulnerabilities.  We  address  and  compare  a  number 
of  frequency  domain  methods  for  suppressing  stationary 
and  slowly  sweeping  narrow-band  interferes.  Frequency- 
domain  interference  excision  techniques  use  the  FFT  to 
identify  and  excise  undesired  spectral  components.  The 
technique  is  known  as  hole-punching,  a  spectral  gating 
operation  that  sets  to  zero  those  spectral  components  ex¬ 
ceeding  specified  threshold  levels. 

The  FFT  of  a  signal  often  exhibits  high  spectral  side-lobes 
due  to  the  finite  aperture  processed  interval.  These  side- 
lobes  represent  spectral  splatter  and  behave  much  like  a 
broad  band  interference  source  by  contributing  energy  to 
all  FFT  spectral  bins.  Suppressing  energy  in  the  main-lobe 
spectral  bins  in  not  adequate  due  to  the  residual  interfer¬ 
ence  distributed  through  the  spectrum  by  these  side-lobes. 
Spectral  splatter  by  high  side-lobes  is  traditionally  called 
spectral  leakage.  We  demonstrate  that  the  spectral  leakage 
that  survives  the  hole  punching  operation  severely  de¬ 
grades  the  performance  of  the  excision  process. 

Excision  performance  is  improved  substantially  with  a  set 
of  algorithms  designed  to  reduce  spectral  side-lobe  levels. 
We  control  the  spectral  side-lobes  of  the  FFT  by  pre¬ 


processing  the  collected  signal  prior  to  the  transform.  In 
the  spectral  analysis  community  this  control  is  accom¬ 
plished  by  the  use  of  data  windows.  Windows  suppress 
spectral  side-lobes  at  the  expense  of  main-lobe  widening 
and  the  rejection  of  data  near  the  boundaries  of  the  data 
collection  interval.  We  compensate  for  this  signal  rejec¬ 
tion  by  overlapping  and  processing  successive  intervals. 
For  windows  with  good  side-lobe  rejection,  the  required 
overlap  is  75%.  This  overlap  seems  to  compound  the 
processing  task  by  requiring  us  to  perform  transforms  four 
times  as  often.  We  avoid  this  increase  in  workload  by 
processing  all  four  overlapped  windows  simultaneously. 
The  resulting  process  is  called  a  poly-phase  filter  bank, 
and  the  combination  of  the  poly-phase  preprocessor  with 
the  FFT  is  termed  a  poly-phase  transform. 

FFT  BASED  TIME  PROCESSING:  The  Fast  Fourier 
Transform  (FFT)  is  a  class  of  algorithms  that  efficiently 
implement  the  Discrete  Fourier  Transform  (DFT).  The 
DFT,  the  discrete  counterpart  to  the  standard  continuous 
domain  Fourier  Transform,  is  traditionally  used  to  form 
the  spectral  description  of  a  sampled  data  time  domain 
signal.  In  a  similar  manner,  the  inverse  DFT  (IDFT)  re¬ 
verses  the  process  to  form  the  sampled  data  time  series 
associated  with  a  given  spectral  description. 

Efficient  FFT  algorithms  are  readily  available  when  the 
size  of  the  transform  is  highly  composite,  that  is,  when  the 
length  N  factors  into  many  small  primes  (or  powers  of 
primes).  The  direct  implementation  of  an  N-point  DFT 
requires  complex  operations.  The  computational 
workload  to  perform  the  same  N-point  transform  as  an 
FFT  depends  to  first  order  on  the  number  of  factors  of  N, 
and  to  second  order,  on  the  values  of  the  factors.  Struc¬ 
tures  designed  to  implement  an  FFT  algorithm  range  from 
variants  of  the  Cooley-Tukey  to  variants  of  the  Good- 
Thomas  (prime  factor)  with  overlays  of  the  Winograd  and 
Rader  Transforms.  In  general,  powers  of  two  or  four,  such 
as  512  or  1024  (2^  or  4^)  point  Cooley-Tukey  variants  are 
the  most  popular  due  to  the  regularity  of  their  butterfly 
and  dragonfly  structures  and  associated  memory  mapping. 
A  close  second  in  popularity  are  the  mixed-radix  Cooley- 
Tukey  routines,  which  support  lengths  such  as  500  or 
1000  (2^x5\  or  2^x5^)  points.  These  are  nearly  as  efficient 
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as  the  single  radix  versions  due  to  the  small  number  of 
sub-structures.  The  third  most  popular  are  the  mixed  radix 
prime  factor  routines  which  support  lengths  such  as  504 
or  1023  (2^x3^x7,  or  3x11x31)  points.  These  are  less  effi¬ 
cient  due  to  the  larger  number  of  sub-structures  or  due  to 
the  larger  valued  primes. 


I  vr;.  |  1  i  vlt 


Figure  1.  Relative  Transform  Time  Versus  Length. 

Figure  1  presents  a  comparison  of  the  relative  time  dura¬ 
tion  required  to  perform  an  FFT  of  different  lengths  of 
approximately  IK.  Here  the  relative  times  are  referenced 
to  the  time  required  to  perform  a  1024  point  transform. 
Denoting  this  reference  as  l-time  unit,  a  1023  point  trans¬ 
form  requires  6.7  time  units,  a  1000  point  transform  re¬ 
quires  3.1  time  units,  and  the  (prime  length)  1021  point 
transform  requires  140.1  time  units.  Efficient  transforms 
are  available  for  lengths  other  than  powers  of  2.  These 
non-traditional  transform  sizes  can  offer  certain  system 
advantages. 

The  ease  with  which  we  can  move  our  description  of  a 
time  signal  from  the  time  domain  to  the  frequency  domain 
and  back  again  has  made  the  FFT  a  core  signal  processing 
function  in  the  digital  signal  processing  (DSP)  commu¬ 
nity.  A  common  signal-processing  task  performed  with 
the  FFT  that  takes  advantage  of  its  computational  effi¬ 
ciency  is  that  of  fast  convolution.  In  this  process,  the  time 
signal  is  first  transformed  to  the  frequency  domain  to  ob¬ 
tain  its  spectral  description.  Here  the  spectrum  is  modified 
using  a  point-by-point  spectral  product  to  form  the  spec¬ 
tral  description  of  the  output  time  series.  This  (modified) 
spectral  description  is  processed  by  an  inverse  FFT  to 
obtain  the  modified  time  series.  This  indirect  process  to 
perform  a  convolution  is  akin  to  performing  an  arithmetic 
multiplication  in  log-space  where  the  operation  (addition) 
is  easier  to  implement.  We  compactly  describe  this  proc¬ 
ess  by  the  statement  that  a  convolution  in  the  time  domain 
is  equivalent  to  a  product  in  the  frequency  domain. 


A  note  of  caution  is  called  for  here.  Processing  consisting 
of  a  transform,  a  spectral  product,  and  an  inverse  trans¬ 
form  in  the  continuous  time  domain  results  in  a  linear 
convolution.  The  same  set  of  operations,  consisting  of  the 
DFT,  a  spectral  product,  and  an  IDFT  results  in  a  circular 
convolution.  When  the  desired  result  of  the  transform 
based  processing  is  a  linear  convolution,  a  mechanism  to 
control  time  domain  wrap-around  or  aliasing  must  be 
imbedded  in  the  block  processing  of  the  FFT.  Standard 
control  techniques  include  zero  extending,  overlap  and 
add,  and  overlap  and  discard  processing. 

The  convolution  process  normally  requires  M  complex 
products  per  output  point  to  implement  an  arbitrary  filter 
of  length  M.  When  M  is  sufficiently  large,  the  transform 
based  fast  convolution  process  offers  a  more  efficient 
implementation.  The  computational  cost  of  performing  a 
convolution  through  the  transform  can  be  estimated  as 
follows:  (N/2)  log2(N)  complex  operations  for  the  input 
transform,  N  complex  operations  for  the  spectral  product, 
and  (N/2)  log2(N)  complex  operations  for  the  output 
transform.  Adding  these  contributions,  we  arrive  at  the 
workload  for  the  transform -based  convolver  as  N 
log2(2N)  complex  operations.  Accounting  for  the  circular 
convolution  properties  of  this  process,  M  points  of  the  N- 
output  points  are  not  useful,  thus  the  number  of  useful 
outputs  for  the  transform-based  process  is  N-M.  The 
workload  per  output  point  (in  units  of  complex  operations 
per  output)  is  the  ratio  (N  log2(2N))/(N-M).  When  the 
transform  length  N  is  large  compared  to  the  filter  length 
M,  this  workload  ratio  is  approximately  log2(2N).  When 
N=1024,  this  workload  is  20.  Thus  we  conclude  that  when 
a  convolver  is  operating  with  filter  lengths  exceeding  ap¬ 
proximately  20  points,  the  transform-based  implementa¬ 
tion  offers  a  significant  computational  advantage.  The 
actual  crossover  point  between  the  two  options  may  vary 
due  to  secondary  considerations  such  as  using  filters  with 
real  or  with  symmetric  weights  as  well  as  length  and  effi¬ 
ciency  of  the  selected  transform.  A  safe  criteria  is  that  the 
crossover  point  is  somewhere  between  16  and  32  points: 
when  the  filter  is  shorter  than  16  point,  do  it  directly, 
when  it  is  greater  than  32  point,  use  the  transform,  and  in 
between,  it’s  probably  a  wash. 

HOLE  POKING  WITH  FFT:  Hole  poking  is  a  simpli¬ 
fied  realization  of  a  circular  convolution  process.  Spectral 
modification  of  the  hole-punch  operation  is  likened  to  an 
ideal  filter  applied  raultiplicatively  in  the  frequency  do¬ 
main.  The  equivalent  filter  exhibits  zero-gain  in  the  stop- 
band  and  unity  gain  elsewhere.  It  almost  sounds  too  good 
to  be  true:  a  perfect  filter!  The  reality  of  course  is  that  it  is 
too  good  to  be  true.  Setting  the  spectral  response  of  an 
FFT  bin  to  zero  does  not  guarantee  that  the  filter  re¬ 
sponse  will  be  zero  or  even  adequate  in  the  neighborhood 
of  that  spectral  zero.  In  fact  the  interval  between  spectral 
zeros  exhibits  the  high  side-lobes  associated  with  sin(x)/x 
which  can  be  as  high  as  0.22  (-13  dB)  relative  to  pass- 
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band  gain  of  unity.  A  typical  spectral  response  of  a  hole- 
punch  filter  defined  on  and  between  the  uniform  grid  of 
an  FFT  is  shown  in  figure  2. 
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Figure  2.  Spectral  Sequence:  Input,  Hole  Punch,  and 
Hole-Punched  Output 


The  top  section  of  this  figure  presents  the  spectrum  of  an 
input  sinusoid  positioned  midway  between  FFT  spectral 
bin  centers.  Due  to  this  positioning,  the  spectral  support 
of  the  input  time  signal  is  not  limited  to  the  immediate 
neighborhood  of  its  peaks.  Gating  the  peaks  with  the  hole- 
punch  filter  (shown  in  the  center  figure)  does  not  affect 
the  input  spectral  side-lobes.  As  seen  in  the  bottom  figure, 
the  peak  amplitude  of  the  surviving  side-lobes  is  ap¬ 
proximately  20%  of  the  original  input  spectral  level. 
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Figure  3.  Input  and  Output  Time  Series  from  Traditional 
Hole-Punching  Process  Corresponding  to  Figure  2. 


Figure  3  presents  the  input  time  signal  that  was  presented 
to  the  processing  represented  by  figure  2.  Also  seen  is  the 
output  signal  obtained  from  that  processing.  Note  that 
while  there  is  significant  reduction  in  the  overall  ampli¬ 
tude  of  the  output  signal  there  are  sections  of  the  time- 
response  that  still  has  significant  amplitude.  In  particular 
the  input  signal  amplitude  is  10  and  sections  of  the  output 
signal  still  has  amplitudes  equal  to  10.  Figure  4  presents 
the  spectral  representation  of  a  hole-punch  process  ap¬ 
plied  to  windowed  input  data.  Note  that  in  response  to  the 
window,  the  input  spectrum  has  reduced  side-lobes  with  a 
slightly  wider  main-lobe.  For  the  same  hole-poke  gating 


shown  in  figure  2,  the  residual  response  of  the  gated 
spectrum  has  peak  amplitude  of  0.7  (as  opposed  to  2.2  in 
figure  2).  The  widening  of  the  main-lobe  suggests  the 
gating  function  should  also  be  widened  and  in  an  actual 
system  we  would  indeed  do  that. 
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Figure  4.  Spectral  Sequence:  Windowed  Input,  Hole 
Punch,  and  Windowed  Hole-Punched  Output 
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Figure  5.  Input  and  Output  Time  Series  from  Windowed 
Hole-Punching  Process  Corresponding  to  Figure  4. 


Figure  5  presents  the  input  time  signal  that  was  presented 
to  the  processing  represented  by  in  figure  4.  Also  seen  is 
the  output  signal  obtained  from  that  processing.  This  fig¬ 
ure  should  be  compared  to  the  equivalent  signal  set  shown 
in  figure  3.  Note  the  significant  reduction  in  the  overall 
amplitude  of  the  output  signal  time-response.  We  must 
also  remember  that  the  window  achieves  this  performance 
by  discarding  input  data  near  the  data  boundaries. 

In  the  previous  paragraph  we  mentioned  that  widening  of 
the  spectral  main-lobe  due  to  the  window  suggests  use  of 
a  wider  hole-punch  gating  function.  Widening  the  hole- 
punch  by  one  index  on  each  side  resulted,  for  this  exam¬ 
ple,  in  a  peak  residual  side-lobe  of  0.04  and  of  time  re¬ 
sponse  of  peak  amplitude  0.1.  A  Kaiser-Bessel  window 
with  -40  dB  side-lobes  was  used  for  this  demonstration. 
Operational  systems  require  additional  side-lobe  rejection 
which  causes  an  increase  in  main-lobe  bandwidth. 
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HOLE  POKING  AS  GPS  OVERLAY:  In  the  last  sec¬ 
tion  we  demonstrated  that  a  windowed  transform  hole¬ 
poking  procedure  offered  significant  performance  advan¬ 
tages  relative  to  the  non-windowed  transform.  We  now 
demonstrate  this  advantage  and  observe  an  undesired  re¬ 
duction  in  processing  gain  of  the  compression  code  due  to 
the  window.  We  will  follow  this  with  a  number  of  over¬ 
lapped  processing  techniques  to  purchase  back  the  loss 
due  to  the  windowing. 
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Figure  6:  Binary  Gold  Code  with  40  dB  Jammer 
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Figure  7:  Output  of  Matched  Filter:  No  Preprocessing 


Figure  6  presents  a  time  signal  consisting  of  a  sinusoidal 
jammer  of  amplitude  100  added  to  binary-phase  Gold 
code  of  amplitude  1  for  a  jammer  to  signal  ratio  of  40  dB. 
Figure  7  presents  the  result  of  a  matched  filter  processing 
the  underlying  binary  phase  Gold  Code.  No  attempt  is 
made  here  to  account  for  analog-to-digital  quantization 
effects  or  receiver  noise.  Note  the  little  bump  at  index 
1023  is  the  desired  compressed  response  from  the 
matched  filter.  The  remainder  of  the  response  is  domi¬ 
nated  by  the  additive  interference. 


Figure  8  presents  a  zoomed  section  of  spectra  obtained  by 
the  FFT  processing  of  the  non-windowed  and  windowed 
data.  Note  the  wider  spectral  width  formed  by  the  non- 
windowed  spectrum.  The  figures  also  present  the  spectral 
gating  interval  to  be  removed  by  the  hole-punch  process¬ 
ing, 


it 

•( 


:•! 


Figure  8.  Spectra  of  Input  Jammer  Signal  plus  Gold  Code 
Obtained  without  and  With  Data  Window 


Figure  9  present  the  outputs  of  the  matched  filter  obtained 
after  the  hole  punching  operations.  The  figures  are  the 
result  of  processing  the  non-windowed  and  windowed 
hole-punch  respectively. 
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Figure  9.  Output  of  M  atched  Filter:  Preprocessed  with 
Non-Windowed  and  Windowed  Hole-Punch  Routines 


Note  the  impressive  rejection  capability  of  the  windowed 
hole-punch  process.  The  penalty  paid  by  the  windowed 
scheme  is  the  reduction  in  processing  gain  due  to  the  win¬ 
dowing  of  the  received  code  as  well  as  the  received  jam¬ 
ming  signal.  The  expected  output  level  from  the  matched 
filler  is  1023,  but  the  windowed  response  is  507,  a  loss 
related  to  the  average  value  of  the  applied  window. 


OVERLAPPED  PROCESSING:  The  transform  used  in 
the  hole-punch  routine  is  a  block  process.  The  desired 
result  of  the  processing  is  a  continuous  flow  of  sampled 
data  representing  a  filtered  version  of  the  collected  data. 
Our  next  task  is  to  convert  the  block  process  to  a  seam¬ 
lessly  merged  set  of  blocks  that  emulate  a  continuous 
flow  process.  Figure  10  is  a  visualization  of  the  time  line 
and  possible  partitions  of  data  into  blocks  to  be  processed 
and  merged  by  the  hole-punch  routine.  The  first  option  is 
to  collect  contiguous  blocks,  process  them  and  reassemble 
their  outputs  into  contiguous  blocks.  This  is  similar  to 
laying  bricks.  We  must  exercise  boundary  control  by  ap¬ 
plying  windows  to  each  data  processing  interval.  If  we 
apply  the  windows  to  the  contiguous  intervals  as  shown  in 
the  second  option  the  data  located  where  the  window  ap¬ 
proaches  zero  is  essentially  discarded  by  the  process.  To 
recover  this  data  we  overlap  the  intervals  to  be  processed 
and  then  assemble  the  final  output  by  summing  the  proc¬ 
essed  overlapped  intervals.  This  process  is  called  merged 
overlapped  windows.  These  techniques  is  similar  to  the 
overlap-and-add  process  used  in  fast  convolution  but  dif¬ 
fers  in  the  use  of  a  data  window  in  each  interval,  a  opera¬ 
tion  that  is  not  used  in  standard  overlap,  and  add  proc¬ 
essing. 
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Figure  10.  Partitioning  the  Time  Line  into  Weighted 
and  Overlapped  Processing  Intervals 


The  intervals  selected  for  processing  can  be  long  or  short, 
relative  to  some  convenient  interval  such  as  the  time  in¬ 
terval  required  for  a  l-K  transform.  The  selection  crite¬ 
rion  is  a  function  of  the  stability  of  the  interference  signal. 
Shorter  processing  intervals  offer  a  small  number  of  wide 
bandwidth  spectral  filters  while  long  intervals  offer  a 
large  number  of  narrow  bandwidth  filters  to  the  hole- 
punch  process.  Figure  11  presents  the  input  and  over¬ 
lapped  windowed  time  interval  partition  presented  to  l-K 
transforms. 


;«  'i  r*  r*  i'*  ^  V:  ‘I  '»  "  *1  ‘‘  *1  •« 

r,-  Jir.  y.  .V.  .V. .  . . 

;.*j  C  X.  ^  Iw."  '.'A* 

iir.i 

-  **  '»  '  **  *r  *»  !*•  •  **  .. 

%  ih  U  fx  ‘  .  •  .  t’  ^  h  f »  c f  \  /  . 

,•  ,  » I  j  <  j  'j  ,  I  i/  i-*  ^  •  n  1 1  ■  ^ 

-  V  'I  i.*  *.*  ‘  ^  " 

ri‘ . . . . . . . -r-— 

C  *A‘-  luC*.  *.*A* 

Figure  11,  Input  and  Windowed  and  Overlap  Partitioned 
Data  for  Successive  l-K  Transforms 


The  hole-punch  algorithm  processes  each  interval  and  the 
resulting  time  series  is  assembled  on  the  time  line  as 
overlapped  contributions  to  the  output  time  series.  These 
overlapped  intervals  are  merged  by  simple  addition  and 
are  presented  to  the  matched  filter  processing  for  subse¬ 
quent  Doppler-removal  and  Gold-code  compression.  Fig¬ 
ure  12  presents  the  real  part  of  the  output  signal  over  a 
section  of  the  time  line  obtained  from  the  overlapped, 
windowed,  and  merged  processing.  The  input  data  here 
was  a  single  interfering  signal  of  amplitude  100  centered 
between  2-FFT  bins.  The  binary  phase  Gold  code  can 
easily  be  seen  on  the  residual  interference  of  approximate 
amplitude  1.  This  represents  a  100-to-l  reduction  in  am¬ 
plitude  to  the  interfering  signal.  This  is  the  expected  level 
of  suppression  that  can  be  obtained  from  the  Hann  win¬ 
dow  selected  for  this  demonstration.  Additional  suppres¬ 
sion  can  be  obtained  with  better  windows  but  we  would 


require  additional  overlap  on  the  order  of  75%  rather  than 
the  50%  used  here. 
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Figure  12.  Time  Response  over  Section  of  Time  Line 
Showing  Spreading  Code  and  Residual  Interference 


Figure  13  presents  the  output  of  the  matched  filter  after 
processing  the  overlapped,  windowed,  and  merged  data 
intervals  with  the  hole-poking  process.  As  expected,  the 
overlapped  processing  recovered  the  processing  gain  loss 
incurred  by  applying  the  window  to  the  data  as  well  as  to 
the  interference  signal.  Note  that  the  peak  output  level  of 
the  filter  is  expected  to  be  the  code  length  of  1023  and  for 
this  example  the  output  was  measured  at  1022. 
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Figure  13.  Matched  Filter  Output  after  Pre-Processing  by 
Windowed,  Overlapped,  and  Merged  Hole-Puncher 


The  next  set  of  figures  demonstrates  the  performance  of 
the  overlapped,  windowed,  and  merged  processing  using 
shorter  transform  intervals.  Figure  14  presents  the  input 
and  overlapped  windowed  time  interval  partitions  pre¬ 
sented  to  256  point  transforms. 
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Figure  14.  Input  and  Windowed  and  Overlap  Partitioned 
Data  for  Successive  0.25-K  Transforms 


These  transforms  are  one-fourth  the  length  of  the  previous 
example  consequently  it  takes  5  overlapped  intervals  to 
span  the  l-K  interval  of  the  Gold  code  as  opposed  to  3 
overlapped  intervals  for  the  previous  example. 

Figure  15  presents  the  real  part  of  the  output  signal  over  a 
section  of  the  time  line  obtained  from  the  window,  over¬ 
lap,  and  merge  processing.  The  input  data  here  was  the 
same  as  that  used  in  the  previous  example  that  differs 
from  this  one  only  in  the  length  of  the  processed  intervals. 
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Here  the  binary  phase  Gold  code  can  easily  be  seen  on  the 
residual  interference  of  approximate  amplitude  2.0.  This 
is  slightly  less  suppression  than  that  obtained  from  the 
larger  transform.  This  is  due  to  the  wider  bandwidths  of 
the  equivalent  filters  presented  to  the  hole-punch  routine. 


that  the  output  of  the  N-point  transform  should  be  sam¬ 
pled  once  for  every  N  input  points  which  is  the  non- 
overlapped  processing  scenario.  The  rectangle  window 
has  unacceptably  high  side-lobes  of  -13  dB  relative  to 
peak  response. 
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Figure  15.  Time  Response  over  Section  of  Time  Line 
Showing  Spreading  Code  and  Residual  Interference 


Figure  16  presents  the  output  of  the  matched  filter  after 
processing  the  shorter  overlapped,  windowed,  and  merged 
data  intervals  with  the  hole-poking  process.  Here  too,  the 
overlapped  processing  purchased  back  the  loss  we  in¬ 
curred  by  applying  the  window  to  the  data  as  well  as  to 
the  interference  signal.  For  this  example  the  output  was 
measured  at  1023.  The  small  increase  in  residual  interfer¬ 
ence  obtained  when  using  the  shorter  transform  intervals 
had  an  insignificant  effect  on  the  side-lode  levels  ob¬ 
served  in  the  output  of  the  matched  filter. 
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Figure  16.  Matched  Filter  Output  after  Pre-Processing  by 
Windowed,  Overlapped,  and  Merged  Hole-Puncher 


WINDOW  CONSIDERATIONS:  We  have  demon¬ 
strated  the  value  of  windowing  to  control  the  spectral 
side-lobe  response  of  the  hole-punching  algorithm.  We 
have  seen  the  need  to  overlap  the  windowed  interval  to 
avoid  rejecting  data  contributions  in  the  time  intervals 
over  which  the  window  amplitude  goes  to  zero.  The  data 
gathered  in  the  overlapped  intervals  is  processed  in  the 
individual  intervals  as  if  they  were  uncorrelated.  Of 
course  the  intervals  are  correlated  and  we  acknowledge 
that  correlation  by  re-merging  the  overlapped  intervals  by 
adding  the  properly  delayed  and  aligned  responses  to  ob¬ 
tain  the  composite  response.  This  overlap  was  described 
for  50%  overlap  in  the  text  associated  with  figure  10  and 
was  demonstrated  in  the  text  associated  with  figures  11 
and  14.  As  mentioned  earlier,  for  good  side-lobe  control 
the  overlap  is  75%.  An  intuitive  argument  justifying  this 
75%  overlap  is  based  on  the  Nyquist  criterion  and  the 
main-lobe  bandwidth  of  the  window  we  use  in  the  hole- 
poke  process.  A  rectangle  window  of  length  N  applied  to 
dau  collected  at  sample  rate  has  a  two-sided  bandwidth 
of  fj/N.  The  Nyquist  rate  for  data  with  this  bandwidth  is 
f,/N  which  is  l/Nth  of  the  input  sample  rate.  This  suggests 


We  control  the  side-lobes  by  applying  a  window  to  the 
data  prior  to  the  N-point  FFT.  If  we  apply  the  Hann  or  the 
Hamming  window  to  realize  -32  dB  or  -44  dB  peak  side- 
lobe  level  we  increase  the  main-lobe  bandwidth  by  a  fac¬ 
tor  of  two  to  2  fs/N  or  fs/(N/2).  The  Nyquist  rate  for  this 
windowed  data  is  f,/(N/2)  which  is  l/(N/2)  of  the  input 
sample  rate.  This  suggests  that  the  output  of  the  N-point 
transform  should  be  sampled  once  for  every  N/2  input 
points  which  is  the  50%  overlap-processing  scenario. 

To  effect  better  side-lobe  control  we  must  select  windows 
with  better  side-lobes.  Good  contenders  are  the  4-term 
Blackman-harris  window  and  the  Kaiser-Bessel  window 
with  time-bandwidth  parameter  B=ll.  These  windows 
have  peak  side-lobes  of  -94  dB  and  -100  dB  respectively. 
They  also  have  a  mainlobe  width  (measured  from  peak  to 
-100  dB  level)  of  4  f^/N  or  f*/(N/4).  Following  the  argu¬ 
ment  we  used  in  the  previous  paragraph,  the  Nyquist  rate 
for  this  windowed  data  is  fs/(N/4).  This  suggests  that  the 
output  of  the  N-point  transform  should  be  sampled  once 
for  every  N/4  input  points  which  is  the  rationale  for  75% 
overlap  processing.  Figure  17  shows  the  partition  and  the 
processing  flow  for  the  75%  overlapped  intervals.  In  par¬ 
ticular  we  see  here  the  time  locations  of  the  partition  and 
merging  processes. 

We  note  that  any  time  interval  of  length  N/4  samples  in 
the  input  data  steam  contributes  to  4-processing  intervals 
and  similarly  any  time  interval  of  length  N/4  samples  in 
the  output  data  stream  is  the  result  of  contributions  from 
4-processing  intervals.  We  also  note  that  the  four  con¬ 
tributors  to  and  from  the  processing  tasks  are  weighted  by 
different  quarter  sections  of  the  window  and  as  each 
fourth  of  the  input  interval  moves  through  the  window  in 
successive  overlapped  intervals  it  is  eventually  weighted 
by  each  of  the  four  fourths.  Specifically,  if  we  examine 
the  banded  section  of  the  overlapped  intervals  shown  in 
figure  17  we  see  the  following.  When  a  new  one-quarter 
input  interval  (of  length  N/4  samples)  is  brought  into  the 
processing  interval  it  is  weighted  by  the  first  fourth  of  the 
window.  When  the  next  quarter  is  introduced  it  is 
weighted  by  the  first  fourth  of  the  window  while  the  pre¬ 
vious  quarter  interval,  previously  the  newest  quarter  and 
subsequently  weighted  by  the  first  fourth  of  the  weighting 
function,  is  now  weighted  by  the  second  fourth.  In  like 
manner  as  each  new  quarter  interval  is  introduced  into  the 
processing  interval  the  previous  intervals  slide  into  later 
sections  of  the  data  window.  We  can  think  of  the  intervals 
as  sliding  through  the  window  in  much  the  same  way  that 
sampled  data  slides  through  a  FIR  filter. 
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Figure  17.  Schematic  Representation  Standard  Processing  of  Input  Interval  Partition,  Windowing,  Processing,  and  Merging 

of  Processed  Data  Intervals  to  Form  Data  Set  in  Final  Output  Series 


An  interesting  awareness  concerning  constraints  on  the 
weighting  function  is  related  to  the  question,  “What  if  the 
processing,  indicated  by  the  processing  blocks  in  figure 
17,  were  simple  identity  operators”?  That  is,  if  we  did 
nothing  to  the  input  signal  except  partition  and  apply  po¬ 
sition  and  partition  dependent  weights  could  we  perfectly 
reconstruct  the  signal?  The  answer  is  yes  if  the  folded 
summation  of  the  four  segments  of  the  weighting  function 
is  a  constant.  A  sufficient  condition  for  perfect  recon¬ 
struction  is  that  the  weighting  term  be  odd  symmetric 
about  the  folding  positions  at  N/4  and  3N/4.  The  symme¬ 
try  condition  just  noted  is  illustrated  in  figure  18. 


Figure  18.  Merging  Sections  of  4-to-l  Overlapped 
Segmented  Weighting  Function  by  Simple  Summation 


The  impact  of  the  symmetry  constraint  just  noted  is  that 
the  spectral  components  of  the  window  are  restricted  to 
DC  and  the  fundamental  of  one  cycle  per  interval.  The 
window  must  be  a  raised  cosine  such  as  the  Hann  or 
Hamming  weights.  If  other  weighting  functions  are  used 
we  will  have  imperfect  reconstruction  of  the  time  series 
which  manifests  itself  as  periodic  modulation  of  the  time- 
domain  envelope.  We  can  elect  to  live  with  the  recon¬ 
struction  artifacts  or  to  suppress  it  with  additional  post¬ 
processing  by  a  second  filter  equivalent  to  a  deterministic 
equalizer. 


POLYPHASE  TRANSFORM:  We  now  take  a  different 
approach  to  the  windowed  transform  processing  required 
to  do  the  spectral  hole-punching.  Till  now,  our  scenario 
has  been,  select  an  interval  and  transform  it.  The  side-lobe 
response  of  the  default  rectangle  window  is  unacceptable 
so  we  apply  a  weighting  function  with  good  spectral  side- 
lobes  to  the  data  interval.  The  use  of  the  window  in¬ 
creased  the  main  lobe  width  of  the  transform,  an  effect 
that  we  turn  to  our  advantage  in  a  later  discussion.  The 
window  also  discarded  data  near  the  boundary  conditions 
of  the  interval.  We  responded  to  this  loss  with  overlapped 
windows  and  a  merging  scheme.  We  now  revisit  the 
comment  that  applying  the  window  increased  the  main- 
lobe  spectral  width.  When  we  apply  a  good  window,  the 
increase  in  spectral  width  is  a  factor  of  4.  Consequently 
the  spectral  measurements  performed  by  the  FFT  are 
highly  correlated,  hence  represent  redundant  data.  Prior 
to  applying  the  window,  the  spacing  between  spectral 
centers  in  the  N-point  FFT  was  f^/N  and  the  bandwidth  of 
the  equivalent  filter  (the  transform  of  the  window)  was 
also  fs/N,  a  good  match!  After  we  apply  the  window  the 
spacing  between  spectral  centers  in  the  N-point  FFT  is 
still  fj/N  but  the  bandwidth  of  the  equivalent  filter  (the 
transform  of  the  new  window)  is  now  4fj:/N,  a  poor 
match!  It  can  be  shown  that  the  summation  of  the  over¬ 
lapped  intervals  represents  additional  processing  that  re¬ 
duces  the  effective  bandwidth  of  these  filters  while  main¬ 
taining  the  desired  good  side-lobe  levels.  We  can  effect 
the  same  results  with  less  computational  effort  by  rede¬ 
fining  the  interval  length  as  well  as  the  window.  Our 
thinking  process  is  we  first  apply  the  window.  The  win¬ 
dow  widens  the  main-lobe  spectral  width  to  4fs/N.  We 
compensate  for  the  increased  numerator  by  offering  a 
matching  increased  denominator  by  building  a  window  of 
length  4N  so  that  the  raainlobe  bandwidth  is  returned  to 
fjj/N  while  realizing  the  desired  good  side-lobe  levels. 

The  first  problem  we  encounter  is  that  the  transform  is 
now  four  times  the  original  size.  We  must  also  remember 
we  still  have  to  use  the  overlapped  intervals  to  avoid  dis¬ 
carding  data.  Let’s  repair  the  first  problem.  If  we  use  a 
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window  of  length  4N  in  a  transform  of  length  4N  we  ob¬ 
tain  spectral  centers  separated  by  f,/4N.  We  originally 
selected  the  processing  interval  of  length  N  because  spec¬ 
tral  spacing  of  f*/N  was  appropriate  to  the  signal  process¬ 
ing  task.  We  also  know  that  the  spectral  measurements  are 
highly  correlated  due  to  widening  of  the  mainlobe  width. 
That  is  the  core  of  our  fix.  We  choose  to  sample  the  out¬ 
put  of  the  transform  not  at  spectral  spacing  of  fs/4N  but 
rather  at  spacing  fs/N  so  that  the  spacing  and  the  filter 
width  once  again  represent  a  good  match.  How  do  we  do 
this  spectral  resampling?  We  simply  skip  samples,  taking 
every  fourth  one  that  the  4N-point  transform  would  offer 
us.  This  is  shown  in  (1)  where  (la)  presents  the  transform 
points  with  poor  side-lobes  but  with  bandwidth  and  spac¬ 
ing  fj/N  obtained  with  an  N-point  FFT.  Equation  (lb) 
presents  the  summation  required  to  obtain  the  same  f,/N 
bandwidth  with  better  side-lobes  due  to  the  4N-point  win¬ 
dow.  Finally  (Ic)  and  (Id)  presents  the  spectral  down- 
sampling  from  fs/4N  to  f^/N  by  computing  every  fourth 
output  point. 
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Noting  that  the  exponential  kernel  in  (Id)  is  periodic  in  N 
samples,  we  can  rewrite  (Id)  to  obtain  (2).  We  recognize 
that  (2)  is  the  sum  of  four  transforms  each  weighted  by 
sections  of  the  window  identified  in  figure  17.  The  data  in 
the  intervals  shown  in  (2)  are  not  the  same  as  that  indi¬ 
cated  in  figure  17.  We  could  perform  four  short  trans¬ 
forms  of  the  segmented  and  weighted  data  and  then  merge 
them  as  suggested  in  figure  17. 
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Since  the  transform  is  a  linear  operator  we  can  rearrange 
the  sum,  and  perform  a  transform  of  a  sum  rather  than  a 
sum  of  transforms.  This  is  shown  in  (3)  where  the  second 
summation  in* (3)  that  computes  P4(n)  is  the  poly-phase 
preprocessor  that  merges  the  weighted  segments  prior  to 
the  transform  as  opposed  to  after. 
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The  impact  of  this  relationship  is  quite  remarkable.  We 
can  perform  non-overlapped  N-point  Transforms  on  75% 
overlapped  and  windowed  sequences  of  length  4N.  This 
is  shown  in  figure  19. 
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Figure  19.  Processing  of  4-to-l  Overlapped  and  Win¬ 
dowed  4N  Point  Intervals  with  Non-Overlapped  N-Point 
Transforms 

The  process  of  wrapping  a  4-N  point  weighted  data  set 
into  an  N-point  Transform  is  called  a  poly-phase  Trans¬ 
form.  The  preprocessing  is  often  called  data  folding  but  it 
is  in  fact  a  data  wrapping  as  can  seen  in  figure  18.  It  is 
intentional  time  domain  aliasing  related  to  the  4-to-l 
downsampling  conducted  in  (1).  Poly-phase  is  a  common 
process  used  in  communication  systems  to  build  high 
performance  spectrum  analyzers  and  orthogonal  fre¬ 
quency  division  multiplexers.  The  process  requires  a 
poly-phase  filter  bank  which  pre-process  short  sequences 
of  data  in  each  bank  and  then  passes  the  processed  data 
vectors  to  the  transform.  This  pre-processor  is  the  second 
equation  in  (3). 
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After  modifying  the  spectrum  by  the  hole-punching  proc¬ 
ess  we  have  to  return  to  the  time  domain.  In  the  previous 
version  of  this  process  we  returned  with  overlapping 
processed  windows  which  had  to  be  merged  by  addition 
of  the  separate  data  sets.  In  this  method  the  data  was 
merged  prior  to  the  transform  so  we  ask  “Are  we  fin¬ 
ished”?  Do  we  simply  take  the  inverse  Transform  and 
merge  the  intervals  by  juxtapositioning  them?  Not  quite! 
We  have  to  perform  the  inverse  of  the  input  procedure.  As 
expected,  the  inverse  entails  the  inverse  transform  but  it 
also  includes  another  poly-phase  filter  bank  that  forms 
weighted  sums  across  successive  intervals  to  invert  the 
folding  process  that  occurred  in  the  input  process.  The 
inversion  would  be  perfect  if  the  window  had  the  sym¬ 
metries  described  earlier.  The  windows  or  filter  weights 
we  use  to  control  side-lobe  levels  contain  additional 
spectral  components  so  that  the  inverse  of  the  input  proc¬ 
ess  does  not  result  in  perfect  time  domain  reconstruction. 
This  can  be  easily  verified  by  skipping  the  forward  and 
inverse  transform  and  simply  have  the  input  poly-phase 
filter  talk  to  the  output  poly-phase  filter. 


The  reason  we  do  not  have  perfect  reconstruction  is  that 
we  have  replaced  the  Dirichlet  kernels  (the  periodic 
sin(x)/x)  equivalent  filters  obtained  by  the  rectangle 
weight  with  better  filters.  The  summation  over  the  fre¬ 
quency  domain  of  the  translated  Dirichlet  kernels  (the 
transform  center  frequencies)  results  in  a  perfectly  flat  or 
constant  spectrum.  The  inverse  transform  of  this  constant 
spectrum  is  the  impulse  response  of  the  system,  which 
happens  to  be  a  single  impulse.  If  we  perform  the  sum¬ 
mation  of  spectral  filters  formed  by  our  windowed  and 
down  sampled  process  we  find  a  periodic  ripple  or  scal¬ 
loping  of  the  spectrum.  This  ripple  in  the  frequency  do¬ 
main  is  observed  in  the  time  domain  as  pre-and  post  ech¬ 
oes  separated  by  the  transform  length.  Thus  an  impulse 
applied  to  the  input  of  the  poly-phase  filter  bank  and  then 
extracted  from  the  output  bank  would  be  accompanied  by 
three  post  echoes  and  three  pre  echoes  (for  4-to-l  over¬ 
lap).  Only  two  of  the  echoes  are  significant  and  they  can 
be  reduced  to  any  desired  level  by  additional  post  proc¬ 
essing  with  separate  deterministic  equalizers  after  each 
poly-phase  output  stage.  What  we  have  accomplished  by 
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Figure  20.  Impulse  Response  of  Composite  Input  and 
Output  Poly-phase  Filter  Banks  and  Equalized  Banks 


this  sequence  of  operations  is  the  replacement  of  the  out¬ 
put  poly-phase  stages  with  a  new  set  obtained  by  con¬ 
volving  each  stage  with  its  deterministic  equalizer.  Figure 
20  presents  the  impulse  response  of  the  standard  poly¬ 
phase  filter  bank  and  the  impulse  response  of  the  equal¬ 
ized  set. 


Figure  21  presents  the  time  series  obtained  from  the 
equalized  poly-phase  filter  bank.  The  input  interference 
was  the  same  sinusoid  of  amplitude  100  used  for  earlier 
windowed  hole-poking  demonstrations.  Note  the  inter¬ 
fering  signal  is  completely  suppressed  and  the  binary 
phase  Gold  code  can  be  seen  in  the  output  series.  There  is 
slight  distortion  in  the  output  series  due  to  the  residual 
system  echoes  remaining  from  imperfect  equalization. 


output  dsoa  from  eq  jaiizBCi  p:](vphase  f  ter  bate 


Figure  21.  Time  Response  over  Section  of  Time  Line 
Showing  Spreading  Code  and  Residual  Interference 


Figure  22  presents  the  output  of  the  matched  filter  after 
processing  the  output  of  the  equalized  poly-phase  filter 
bank.  As  expected  the  full  compressed  signal  is  present 
along  with  the  side-lobe  levels  of  the  despread  code.  Es¬ 
sentially  no  degradation  results  in  the  expected  signal 
structure  due  to  the  hole-poking  and  the  overlapped  poly¬ 
phase  processing.  The  expected  peak  output  level  is  1023, 
the  measured  response  was  1004. 
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Equalized  Poly-Phase  Transform  Based  Hole-Puncher 


Finally,  for  completeness,  figure  23  presents  the  output  of 
the  matched  filter  after  processing  the  output  of  the  non- 
equalized  poly-phase  filter  bank.  The  compressed  signal 
is  present  along  with  the  side-lobe  levels  of  the  despread 
code.  We  also  see  degradation  of  the  expected  signal  due 
to  the  system  echoes.  In  particular  note  the  two  echoes  at 
time  indices  4613  and  5125.  The  Main  peak  response  is  at 
location  4869,  so  echoes  are  where  they  were  expected  at 
±  N  (transform  length  N=256).  The  expected  peak  output 
level  is  1023,  the  measured  response  was  only  868. 
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In  this  paper  we  presented  one  powerful  application  of  the 
transform  to  excise  strong  spot  interference.  The  trans¬ 
form  can  efficiently  perform  convolutions  and  correla¬ 
tions.  In  our  proposed  next  generation  receiver  architec¬ 
ture  we  have  access  to  the  transform  which  we  use  to  per¬ 
form  several  tasks  that  are  not  implemented  in  current 
receivers.  The  transform  assists  in  carrier  acquisition  by 
defining  Doppler  lines.  It  also  assists  in  code  acquisition 
by  performing  the  entire  code  phase  hypothesis  test  in  a 
single  look  at  the  data.  The  transform  is  also  used  to  in¬ 
vert  multi-path  channels  by  absorbing  the  channel  inverse 
in  the  transform  of  the  spreading  code’s  matched  filter. 
This  work  will  be  documented  in  forthcoming  publica¬ 
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